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Abstract 



We calculate the thermodynamic functions of a hot gluon plasma to leading 
order in hard-thermal-loop (HTL) perturbation theory. Effects associated 
with screening, gluon quasiparticles, and Landau damping are resummed to 
all orders. The ultraviolet divergences generated by the HTL propagator 
corrections can be canceled by a counterterm that depends on the thermal 
gluon mass parameter. The HTL thermodynamic functions are compared to 
those from lattice gauge theory calculations and from quasiparticle models. 
For reasonable values of the HTL parameters, the deviations from lattice 
results for T > 2Tc have the correct sign and roughly the correct magnitude 
to be accounted for by next-to-leading order corrections in HTL perturbation 
theory. 
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I. INTRODUCTION 



Relativistic heavy-ion collisions will soon allow the experimental study of hadronic mat- 
ter at energy densities that will probably exceed that required to create a quark-gluon 
plasma. A quantitative understanding of the properties of a quark-gluon plasma is essen- 
tial in order to determine whether it has been created. Because QCD, the gauge theory 
that describes strong interactions, is asymptotically free, its running coupling constant as 
becomes weak at sufficiently high temperatures. This would seem to make the task of under- 
standing the high-temperature limit of hadronic matter relatively straightforward, because 
the problem can be attacked using perturbative methods. Unfortunately, a straightforward 
perturbative expansion in powers of as does not seem to be of any quantitative use even 
at temperatures that are orders of magnitude higher than those achievable in heavy-ion 
collisions. 

The problem is evident in the free energy of the quark-gluon plasma, whose weak- 
coupling expansion has been calculated through order a^/^ HhSl- An optimist might hope to 
use perturbative methods at temperatures as low as 0.3 GeV, because the running coupling 
constant as(2vrT) at the scale of the lowest Matsubara frequency is about 1/3. Unfortunately, 
the weak-coupling expansion seems to diverge badly even at much higher temperatures. For 
a pure-glue plasma, the first few terms in the weak-coupling expansion are 
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where J-'ideai = — (87r^/45)T^ is the free energy of an ideal gas of massless gluons and ag = 
as(/i4) is the running coupling constant in the MS scheme. In Fig. |l|, the free energy is shown 
as a function of T/Tc, where is the critical temperature for the deconfinement transition. 
The weak-coupling expansions through orders a^, a^^^, a^, and a^^^ are shown as bands that 
correspond to varying the renormalization scale /i4 by a factor of two from the central value 
fi4 = 2'kT. As successive terms in the weak-couphng expansion are added, the predictions 
fiuctuate wildly and the sensitivity to the renormalization scale grows. Of course, because 
of asymptotic freedom, the first few terms in the weak-coupling expansion will appear to 
converge at sufficiently high temperature. However, this occurs only at temperatures orders 
of magnitude larger than T^. For example, the a^^^ term is smaller than the as term only if 
as is less than about 1/20, which corresponds to a temperature greater than about lO^Tc. It 
is clear that a reorganization of the perturbative series is essential if perturbative calculations 
are to be of any quantitative use at temperatures accessible in heavy-ion collisions. 

The poor convergence of the perturbative series is puzzling, because lattice gauge theory 
calculations indicate that the free energy JF of the quark-gluon plasma can be approximated 
by that of an ideal gas unless the temperature T is very close to Tc @J^. In Fig. [l|, we also 
show the lattice results for the free energy of pure-glue QCD from Boyd et al. |Q. The free 
energy approaches that of an ideal gas of massless gluons from below as the temperature is 
increased and the deviation is less than about 25% if T is greater than 2Tc. 
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FIG. 1. The free energy for pure-glue QCD as a function of T/Tc- The weak-coupHng expan- 
sions through orders a^, a^^^, a^, and Oi'J'^ are shown as bands that correspond to varying the 
renormahzation scale /i4 by a factor of two. The lattice results from Boyd et al. Q are shown as 
diamonds. 



There are many possible ways to reorganize the perturbation series for the free energy in 
order to improve its convergence. One possibility is to apply Pade approximation methods 
to the expansion in a^"^ Q]. This does improve the convergence somewhat, but it is at best 
a recipe with little physical motivation. There are also technical problems associated with 
the fact that the weak-coupling expansion is not simply a power series in aj/^. Effective 
field theory methods can be used to show that it has the general form 

oo 

^QCD = T4Ecn(loga.X/2, (2) 

n=0 

where the coefficients c„(logas) are polynomials in loga<j 0. The first few coefficients 
can be calculated using perturbative methods and are given in (|l]). Beginning at order 
a^, nonperturbative methods are required to calculate some of the coefficients. Another 
problem with the Pade approximation method is that it is applicable only if several terms in 
the perturbative series are known. This essentially limits its applicability to the free energy 
JF and to the other thermodynamic functions that can be obtained from JF by differentiation. 
Almost all calculations involving signatures of the quark-gluon plasma have been carried out 
only to leading order. The only exception is the production of hard dileptons which has been 
calculated to next-to- leading order in [Q. 

Another approach to the problem is based on the observation that the large corrections of 
order a^J'^ and a^J'^ in (|lD arise from the momentum scale gT. Effective field-theory methods 



3 



can be used to separate the momentum scale T from the lower scales gT and g^T. The 
resulting effective field theory, dimensionally- reduced QCD (DRQCD), is an SU{3) gauge 
theory with adjoint scalars in three euclidean dimensions. The construction of DRQCD can 
be used to separate the free energy into contributions from the scale T and from lower scales: 

^(T) = J-o(T, as) + T/i(m|, gl, ...). (3) 

The function JFq and the parameters m|;, g'^, ... of the effective field theory can be computed 
as power series in as{27rT) by matching the effective field theory with thermal QCD. For 
pure-glue QCD, the parameters of DRQCD at leading order in as are m|; = g^T"^, g% = g^T. 
The function fi can be expanded in powers of gl^/iriE and the coefficients of the first three 
terms are known analytically P]. The first term gives the large a^J"^ correction in (|I]). The 
third term contributes — 830(as/7r)^/^ to the expansion in (|I]) and therefore accounts for most 
of the large a^J'^ correction. Thus the poor convergence properties of the QCD perturbative 
series seems to arise from a breakdown of perturbation theory in DRQCD. This suggests 
that nonperturbative methods should be used to compute the second term in as a 
function of the parameters m|;, g\., ... of DRQCD. The only systematic nonperturbative 
method currently available is Monte Carlo simulations of lattice DRQCD. This method has 
been developed by Kajantie, Rummukainen, and Shaposhnikov 0, and used to compute 
the Debye screening mass for thermal QCD It has not yet been applied to the free 



energy. One disadvantage of this method is that it relies in an essential way on the analytic 
continuation of the quantum field theory to imaginary time. It therefore cannot be used to 
calculate signatures of the quark-gluon plasma that involve real-time processes. 

The poor convergence of the weak-coupling expansion for the free energy is not specific 
to QCD. The free energy of a massless scalar field theory with a 0^ interaction has also 
been computed to order g^ O]. The series seems to converge only if the coupling constant 
g is extremely small. The largest corrections, which include the g^^"^ term, come from the 
momentum scale gT. This suggests that the source of the convergence problem could be 
similar for QCD and the massless scalar theory. 

There have been several other attempts to improve the convergence of the free energy for 
the massless scalar field theory One of the most successful approaches is "screened 



perturbation theory" developed by Karsch, Patkos, and Petreczky [|T2[. This approach 



can be made more systematic by using the framework of "optimized perturbation theory" , 
which has been applied by Chiku and Hatsuda to a scalar field theory with spontaneous 
symmetry breaking . A local mass term proportional to 0^ is added and subtracted from 
the lagrangian, with the added term included nonperturbatively and the subtracted term 
treated as a perturbation. The renormalizability of the 0^ term guarantees that all ultraviolet 
divergences generated by the mass term can be systematically removed by renormalization. 
When the free energy is calculated using screened perturbation theory, the convergence of 
successive approximations to the free energy is dramatically improved. 

Conventional perturbation theory is essentially an expansion about an ideal gas of mass- 
less particles. Screened perturbation theory is a reorganization of perturbation theory such 
that the expansion is about an ideal gas of quasiparticles with a temperature-dependent 
mass. Empirical evidence that such a reorganization of perturbation theory might be useful 
for QCD is provided by the success of quasiparticle models for the thermodynamics of QCD 
at temperatures above T^.. There were several early attempts to fit the thermodynamic 
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functions calculated by lattice gauge theory to those of an ideal gas of massive quarks and 
gluons 0,11^]. In recent years, there have been several significant improvements in the lat- 
tice gauge theory calculations. The thermodynamic functions for pure-glue QCD have been 
calculated with very high precision by Boyd et al. |Q. There have also been calculations with 
Nf = 2 and 4 flavors of dynamical quarks . Motivated by these developments, more quan- 



titative comparisons with quasiparticle models have been carried out by Peshier et al. ||T8 
and by Levai and Heinz These analyses indicate that the lattice results can be fit 

surprisingly well by an ideal gas of massive quarks and gluons with temperature-dependent 
masses mq{T) and mg{T) that grow approximately linearly with T. 

A straightforward application of screened perturbation theory to gauge theories like 
QCD would be doomed to failure, because a local mass term for gluons is not gauge invari- 
ant. There is a way to incorporate plasma effects, including quasiparticle masses, screening 
of the gauge interaction and Landau damping, into perturbative calculations while main- 
taining gauge invariance and that is by using hard-thermal-loop (HTL) perturbation theory. 
Hard-thermal-loop perturbation theory was originally developed to sum up all higher loop 
corrections that are leading order in g for amplitudes having soft external lines with mo- 
menta of order gT pO[. If it is applied to amplitudes with hard external lines with momenta 
of order T, it selectively resums corrections that are higher order in g. This resumma- 
tion is a generalization of screened perturbation theory that respects gauge invariance. It 
corresponds essentially to expanding around an ideal gas of quark and gluon quasiparticles. 

HTL perturbation theory is defined by adding and subtracting hard-thermal-loop cor- 
rection terms to the action pO|- The HTL correction terms are nonlocal, and the resulting 
effective propagators and vertices are complicated functions of the energies and momenta. 
Calculations in HTL perturbation theory are therefore much more difficult than in screened 
perturbation theory. The nonlocality of the HTL correction terms also introduces concep- 
tual problems associated with renormalization. Since these terms are not renormalizable in 
the standard sense, the general structure of the ultraviolet divergences that they generate 
is not known. 

If HTL perturbation theory proves to be tractable and if issues associated with renor- 
malization can be resolved, it could have very important applications. It is a reorganization 
of QCD perturbation theory around a starting point that is essentially an ideal gas of mas- 
sive quasiparticles. In contrast to quasiparticle models, the corrections due to interactions 
between quasiparticles can be calculated systematically, at least through next-to-next-to- 
leading order in HTL perturbation theory. Beyond that order perturbative calculations 
must be supplemented by a nonperturbative method to deal with the magnetic mass prob- 
lem. A significant advantage of HTL perturbation theory over the approach using lattice 
gauge theory calculations in DRQCD is that it can be readily applied to the real-time pro- 
cesses that are the most promising signatures of a quark-gluon plasma. 

In this paper, we calculate the free energy of a hot gluon plasma explicitly to lead- 
ing order in hard-thermal-loop perturbation theory. In spite of the complexity of the HTL 
propagators, their analytic properties can be used to make explicit calculations possible. 
Although complicated ultraviolet divergences arise in the calculation, the most severe diver- 
gences cancel between quasiparticle and Landau-damping contributions and between trans- 
verse and longitudinal contributions. The remaining divergences arise from integrating over 
large three-momentum and can be removed by a counterterm proportional to at the 
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expense of introducing a renormalization scale. With reasonable choices of the renormaliza- 
tion scales, our leading order result for the HTL free energy lies below results from lattice 
QCD for T > 2Tc. However, the deviation from lattice QCD results has the correct sign 
and roughly the correct magnitude to be accounted for by next-to-leading order corrections 
in HTL perturbation theory. 



II. HTL RESUMMATION 

We begin this section by defining screened perturbation theory for a scalar field theory. 
We then define hard-thermal-loop perturbation theory, which is a generalization of screened 
perturbation theory that respects gauge invariance. Finally, we write down a formal expres- 
sion for the free energy at leading order in HTL perturbation theory. 



A. Screened Perturbation Theory 

The lagrangian density for a massless scalar field with a 0^ interaction is 

C = ^d,<pdy-^g'<P' + AC{g'), (4) 

where g is the coupling constant and AC includes counterterms. The conventional pertur- 
bative expansion in powers of g^ generates ultraviolet divergences, and the counterterm AjC 
must be adjusted to cancel the divergences order by order in g'^. At nonzero temperature, 
the conventional perturbative expansion also generates infrared divergences. They can be 
removed by resumming the higher order diagrams that generate a thermal mass of order gT 
for the scalar particle. This resummation changes the perturbative series from an expansion 

1 /2 

in powers of g'^ to an expansion in powers of {g'^) = g. Unfortunately, the perturbative 
expansion for the free energy has large coefficients and appears to be convergent only for 
tiny values of g |[Tl[ . 



Screened perturbation theory, which was introduced by Karsch, Patkos and Pe- 
treczky [|12| , is simply a reorganization of the perturbation series for thermal field theory. It 



can be made more systematic by using a framework called "optimized perturbation theory" 
that Chiku and Hatsuda ||T^ have applied to a spontaneously broken scalar field theory. 



The Lagrangian density is written as 

C = C-^{m^ - + ACs{g\ - W), (5) 

where 8 = \ will be used as a formal expansion parameter and AC^ includes the additional 
counterterms that are required to remove ultraviolet divergences for 6 If we simply set 
(5 = 1, the additional terms in vanish. Screened perturbation theory is defined by taking 
VP? to be of order g^ and 5 to be of order g^^ expanding systematically in powers of g and 
setting 5 = 1 at the end of the calculation. This defines a reorganization of perturbation 
theory in which the expansion is around the free field theory defined by 

Aree = ^5^09^ - \m^\ (6) 
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The effects of the term in (P) are included to all orders, but they are systematically 
subtracted out at higher orders in perturbation theory by the 5m? term in (^). 

This reorganization of perturbation theory generates new ultraviolet divergences, but 
they can be canceled by the additional counterterms in A£s. The renormalizability of 
the lagrangian in guarantees that the only counterterms that are required are df^(j)d'^(j), 
0^, and (p^. At nonzero temperature, screened perturbation theory does not generate any 
infrared divergences, because the mass parameter m? in the free lagrangian (|^) provides an 
infrared cutoff. The resulting perturbative expansion is therefore a power series in g"^ whose 
coefficients depend on the mass parameter m^. 

The parameter in screened perturbation theory is completely arbitrary. To complete 
a calculation in screened perturbation theory, it is necessary to specify as a function 
of and T. Karsch, Patkos, and Petreczky used the solution of a gap equation as their 
prescription for ir?{T). The resulting loop expansion for the free energy appeared to be 
convergent. 

In the weak coupling limit (7^0, the solution to the gap equation for im?{T) approaches 
5f^r^/24. If we use this value for m? and then reexpand the perturbative series for the free 
energy in powers of g, we recover the conventional perturbative expansion in powers of g, 
with its lack of convergence. It is therefore essential to keep the full dependence of m? at 
every order in g"^ and not expand in powers of m/T . 

B. HTL Perturbation Theory 

Renormalized perturbation theory for pure-glue QCD can be defined by expressing the 
QCD lagrangian density in the form 

-^QCD ~ ^'2^'^ {G^uG^'^) + Cgf + £ghost + '^'^QCD, (7) 

where G^i, = d^Ay — d^A^ — ig[A^, A^] is the gluon field strength and A^ is the gluon field 
expressed as a matrix in the SU (3) algebra. The ghost term £ghost depends on the choice of 
the gauge-fixing term £gf. The perturbative expansion in powers of g generates ultraviolet 
divergences, and the counterterms in A£qcd are adjusted to cancel those divergences order 
by order in g. 

Hard-thermal-loop (HTL) perturbation theory is simply a reorganization of the pertur- 
bation series for thermal QCD. The lagrangian density is written as 

>Cqcd = -^^QCD + -^^HTL + A£htl(5', m^- 5ml). (8) 

The HTL improvement term is 

W = -\H- K)Tr (g,. (^^)^ , (9) 

where is the covariant derivative in the adjoint representation, = (1, fi) is a light-like 
four- vector, (. . .)n represents the average over the directions of fi, and 5 = 1 is a formal 
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expansion parameter. The term (|^) is the effective lagrangian that would be induced by 
a rotationally invariant ensemble of colored sources with infinitely high momentum. The 
parameter can be identified with the plasma frequency, or equivalently, with the rest 
mass of a gluon quasiparticle. 

If we set 5 = 1, the coefficient of the HTL improvement term (^) vanishes. HTL 
perturbation theory is defined by taking to be of order and 6 to be formally of order 
g"^, expanding systematically in powers of g, and then setting 5 = 1 at the end of the 
calculation. This defines a reorganization of the perturbative series in which some of the 
effects of the term in are included to all orders but then systematically subtracted 
out at higher orders in perturbation theory by the Sm^ term. If the perturbative corrections 
from the term could be summed to all orders, there would be no dependence on rrig. 
However, any truncation of the expansion in 6 produces results that depend on rrig. 

The term A£htl in ® includes the counterterms required to cancel the new ultra- 
violet divergences generated by the reorganization of the perturbative series. Unlike the 
counterterms A£qcd in ® which are constrained by locality and renormalizability, the 
general structure of the terms in A£htl is unknown. We anticipate that if we use a scale- 
invariant regularization method such as dimensional regularization, the terms in A£htl will 
be polynomials in the parameters g and — Sm^. 

The free lagrangian that serves as a starting point for HTL perturbation theory is 
obtained by setting g = and (5 = in (H). Choosing a covariant gauge-fixing term with 
gauge parameter ^, the free lagrangian is 

1„ . .2 



-|m^lV (^{d,A„ - a„A„) (^^^) {^A, - a.A"')]^ . (10) 

The HTL improvement term generates a self-energy tensor for the gluon that is diagonal in 
the color indices and has the form 

where = (A;o,k) is the Minkowski four-momentum and = — k^. The HTL self- 
energy tensor satisfies the Ward identity K^Jl^^'^ = 0. Because of this Ward identity and the 
rotational symmetry around the axis k, the tensor H^'^ can be expressed in terms of two 
independent functions H^ and H^: 



g^^S'^ - k'y)g^\ (12) 



The longitudinal and transverse self-energies are 



Hi(K) = H°0(ir), (13) 
1 
2 



I[t{K) = = \ (5'^ - Uy) W^{K). (14) 
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The inverse propagator for the gluon is diagonal in the color indices and has the form 



A-^K)"^ = -K^g^'' + - K^K^ + W{K). (15) 

If there are d spatial dimensions, the matrix A~^{K)'^^ has d + 1 eigenvalues: —K'^{k'^ — 
IlL)/k'^, —K'^/^i and the (d — l)-fold degenerate eigenvalue K^ — Ut- The inverse propagator 
for the ghost in a general covariant gauge is Agj^^^^^^ = K'^. 



In the general Coulomb gauge, the gauge- fixing term in ([T0|) is replaced by 
-(1/OTr ((V ■ A)^). The term -{l/^)K^'K'' in the inverse gluon propagator (p!5| ) is re- 
placed by —(l/^)g^^k'^k^g^'^. The d + 1 eigenvalues of A^^(K)^^ are the {d — l)-fold degen- 
erate eigenvalue — 11^ and two other eigenvalues whose product is —k'^ikP' — Yii)/^. The 
inverse propagator for the ghost in this gauge is fc^. The limit ^ — is the strict Coulomb 
gauge defined by the constraint V ■ A = 0. In this gauge, the only propagating modes of 
the gluon are and the d — 1 transverse components of A. The propagators are k^ — Hl 
for Aq and — Ht for the transverse components of A. 



C. HTL Free Energy 

In the imaginary-time formalism, the renormalized one-loop free energy can be written 



as 



^HTL = {Nl - I] 



\f^Ti\ogA-\K) -fjogA-J,,,,{K) + AT 



(16) 



where and Agj^^^^. are the euclidean inverse propagators for gluons and ghosts, and K 
is now a euclidean four-momentum: = (ci;„,k), and = k"^ + u'^. The sum-integral 
in (|l^) represents a dimensionally regularized integral over the momentum k and a sum over 
the Matsubara frequencies u;„ = 27rnT: 

The factor of where is a renormalization scale, ensures that the regularized free 

energy has correct dimensions even for d ^ ?>. The counterterm AjF in (|TBp can be used to 
cancel ultraviolet divergences that have the form of an additive constant in the free energy. 
In the gluon term in (|TB|), the integrand is the sum of the logarithms of the d + 1 



eigenvalues of A~^. After taking into account the cancellations from the ghost term and 
dropping terms that vanish in dimensional regularization, the expression (0) for the HTL 
free energy reduces to 

^HTL = {Nl - 1) [{d - 1)Tt + Tl + A^] , (18) 

where 
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(19) 
(20) 



An identical expression for the free energy is obtained in the general Coulomb gauge. The 
free energies J-'t and J-'l have simple interpretations in the strict Coulomb gauge. They are 
simply the contributions to the free energy from a transverse component of A and from Aq, 
respectively. We therefore refer to them as the transverse and longitudinal free energies. 

At first sight, the calculation of the HTL free energy appears to be a rather daunting 
mathematical problem. The HTL self-energy functions in ( |T9| ) and (pO|) are 



2^ F 



1 + — — log 



2iUr,k 



iuj„ — k 



nL(u;„, k) = 3ml 



2k iujr, — k 



(21) 
(22) 



In order to calculate the HTL free energy, we must evaluate the sum- integrals in ( 19 ) and (^OD 
with a regularization that cuts off the ultraviolet divergences. In the next section, we show 
that this problem can be made tractable by using the analytic properties of the HTL self- 
energies. 



III. HTL FREE ENERGY 



In this section, we reduce the HTL free energy to integrals that can be evaluated numer- 
ically. We first compute the free energy of a free massless boson as a simple example. We 
then compute the transverse and longitudinal free energies, separating them into quasipar- 
ticle and Landau-damping terms and isolating the ultraviolet divergences into integrals that 
can be computed analytically. The most severe divergences cancel, leaving a logarithmic 
divergence that must be canceled explicitly by a counterterm. 



A. Free Massless Boson 



As a warm-up exercise for computing the transverse free energy, we compute the free 
energy of a free massless boson: 

-^boson = ^^^log(ir'). (23) 

In dimensional regularization, the n = term in the sum-integral vanishes. For the same 
reason, we can subtract log/c^ from the integrand for each of the n 7^ terms. The free 
energy becomes 

^boson=^TEXlog^-^. (24) 
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We have introduced a compact notation for the dimensionally regularized integral over mo- 
mentum together with the renormalization scale factor: 



/i 



3-d 



d' 



(25) 



The sum over Matsubara frequencies in ( |24D can be expressed as a contour integral: 

1 



>bo.on -2Jc2mk k^ eP^ 



where the contour C encloses the points ui = iun, 7^ 0, as shown in Fig. 



(26) 



FIG. 2. The free energy of a massless boson can be expressed as an integral over a contour C 



that can be deformed into two contours that wrap around the branch cuts of Iog(A; 



■ UJ 



We can choose the branch cuts of the logarithm to run along the real cj-axis from —00 
to —k and from +k to +00. It is convenient to insert a convergence factor e'''^ {r] 0"*") into 
the integrand. This allows the contour C in Fig. |^ to be deformed into two contours that 
wrap around the branch cuts. Collapsing the contour onto the two branch cuts, making the 
change of variables u —uo in the contribution from the negative branch cut, and taking 
the limit 77 — > whenever it is not needed for convergence, the free energy reduces to 



boson 



+ 



(27) 



Integrating over uj and taking the limit ^ 0, we obtain 

^boson= / Tlog(l-e-'^^)-^ + iA; 



(28) 
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The integral over the last two terms are zero in dimensional regularization. Setting d = 3 
in the integral of the first term in (|28|) and evaluating it analytically, we get 

•^boson ~ rin'^ ' (^'^) 

yu 

Summing over the — 1 color degrees of freedom and the two spin degrees of freedom of 
a transverse gluon, we obtain the free energy of an ideal gas of massless gluons: 

^ideal = - {N^ - l) (30) 



B. Transverse Free Energy 



We now proceed to compute the transverse free energy in ([T9|) . Following the steps that 



led from ( p3D to (p6D, we obtain the contour integral 



T 



duo 

2 Jc 2m Jk 



log- 



+ IlT{—iuJ,k) e 



1' 



(31) 



where e^"^ (?] — > 0"*") is a convergence factor and the contour C encloses the points uj = 
iuinifi 7^ 0. The integrand has logarithmic branch cuts in uj that run along the real cu-axis 
from — oo to —ujT{k) and from +uiT{k) to +cxd, where Uxik) is the quasiparticle dispersion 

.2 



relation for transverse gluons. This dispersion relation satisfies k — uj^ + IiT{—i^T} k) 
or 



0, 



(^2 _ l2 , ^^2^^ 

^ ~ 2 ^ A;2 



T 



2ujrk 



k"^ . ut + k 
log 



LOT 



(32) 



The integrand also has a branch cut in uo running from —k to +k due to the function TIt- 
This branch cut represents the effects of Landau damping. The contour C can be deformed 
into a contour that wraps around the three branch cuts as in Fig. |^. 

We identify the contribution from the branch cuts ending at -^uoxik) as the quasiparticle 
part of jFj. and denote it by J^T,qp- Following the same steps that led to (|2^) , we obtain 



T,qp 



Tlog(l - e 



-Plot' 



1 

H — uot 
2 ^. 



(33) 



We have omitted the term 1/(2,7]) in the integrand of ( P3| ) because its integral vanishes in 
dimensional regularization. We identify the contribution from the contour wrapping around 
the branch cut running from —k to +k as the Landau-damping part of J^t and denote it 
by J-'tm- Collapsing the contour onto the branch cut and expressing it as an integral over 
positive values of a;, we obtain 



T,Ld 



duo 



n Jk Jo 



1' 

T+2, 



(34) 



The angle (pT is 
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FIG. 3. The transverse free energy can be expressed as an integral over a contour C that wraps 



around three branch cuts of log(/c 



+ nr). 



Pt{'^i k) = arctan 



4 



1 -^ T 



(35) 



where = k"^ — uo'^ and L = log [(A; + Ci;)/(A; — u)]. The angle 4>t{^^, k) vanishes as ^ 
and as A; — s> oo and it also vanishes dX uo = k. The complete transverse energy Tt is the 
sum of the quasiparticle term (|33|) and the Landau-damping term (|3^) . The integral of the 
second term in ( P5| ) and the integral of both terms in ( P^ are ultraviolet divergent, but the 
divergences are regularized by the li-dimensional integral over k. In order to extract the 
divergences analytically, we make subtractions in the integrands that render the integrals 
finite in c? = 3 dimensions and then extract the poles in d — 3 from the subtracted integrals. 

We first consider the transverse quasiparticle term (0). The integral of the second term 
is divergent, because the asymptotic form of the transverse dispersion relation for large k 

is m 



, ^ ^rnl 9 / gk^ 
^+4X-32l^r^°^« 



(36) 



Our subtraction should make the integral ultraviolet convergent for d = 3, and it should not 
introduce any infrared divergences. Our choice for the subtracted integral is 

Kk' + Imp - 



(sub) 



T,qp 



1 

2 Jk 



9m^ 



+ 16(^2 + 1^2)3/2 




3m g 



After subtracting this from (|3^), we can take the limit d —>■ 3: 



4vr2 



k'dk 



ut - \ k^ + iml + 



9m^ 



16(A;2 + |m2)3/2 



log 



l{k'' + f m^) 
3m?g 



(37) 



(38) 
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If we impose a momentum cutoff k < A, our subtraction integral (^) has power divergences 
proportional to A*^ and m^A^ and logarithmic divergences proportional to m^logA and 
m^log^A. The quartic divergence is canceled by the usual renormalization of the vacuum 
energy density at zero temperature. Dimensional regular izat ion throws away the power 
divergences and replaces the logarithmic divergences by poles in d — 3. In the limit d ^ 3, 
the individual integrals in (^) are given by ( |A.1| )-( [X7^ ) in the appendix. The result is 



(sub) 



T,qp 



9 4 



ft. 



(2vr) 



1 41og2 
^ + ^ 



31 , 2 , TT^ 1 

+ 21og'2-31og2 + — -- 

e 6 4 



(39) 



where e = {3 — d)/2 and Qd = 27r'^/^/r((i/2) is the angular integral in d dimensions. 

We next consider the Landau-damping term (0). Again we must choose a subtraction 
term that removes the ultraviolet divergences without introducing infrared divergences. Our 
choice is 



(sub) 



T,Ld 



1 



TT Jk JO 



duj 



37rm^ci; 



4^3 g/3a; 



+ 



3'KmlujK^ 



2\2 



|m2)2 



After subtracting this from (|3^) , we can take the limit d —>■ 3: 



2k 



L-1 



(40) 



/J-- ^(sub) 



27r3 



duj- 



,f3LU 



Pdk 



3Tirn?gUJ 



47^3 



doj / k'^dk 



3'KmluK^ 

3 

2 ■■"9 



+ 



4A;3 



Ak^{K^ + lml) ' Sk^{K^ + lmlf \2k 



. (41) 



In both integrals, the first subtraction makes them convergent as /c — ^ cxd with uj fixed. The 
second subtraction in the second integral makes it convergent as A; — ^ oo with uj ^ k. If we 
impose ultraviolet cutoffs a; < A and /c < A on the energy and momentum, the subtraction 
integral (^O]) has a power divergence proportional to m^A^ and logarithmic divergences pro- 
portional to m^T^ log A, m^logA, and m^log^A. The divergence proportional to m^log^A 
cancels against the corresponding divergence in the quasiparticle subtraction integral (p7[). 
The cancellation can be traced to the fact that lixi'^ik) is analytic in the variable uj at 
uj = oo. Dimensional regularization replaces the logarithmic divergences in (^) by poles in 
d — 3 and sets the power divergences to zero. The integrals in (^) are evaluated in the limit 
(i — > 3 in the appendix and given by ( [A.4| )-( [70D and ( |A.9| ). The result is 



^(sub) 
r,Ld 



16 ^ V/i2 



fin 



(2vr) 



1 + 2^ 
.e+'c(-l) 



21og(27r) 



(5) 


' fid 




{2'kY 



2 1og2-l 1 2, 2 22, 10 

+ - bg^2 bg2 + — + — 

3 e 3^ 9^ 6 9 



(42) 



Our final result for the transverse self-energy is the sum of 
that the divergence proportional to l/e^ cancels between ( ^9]) and (|4^) 



Tl) and 



Note 
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C. Longitudinal Free Energy 



We now compute the longitudinal free energy given in (|20|) . We first separate out the 
n = term: 



rL = \T I log ie + 3m2) + \tY.( log - IlL{uJn. k 



(43) 



These two terms can be expressed as a single integral over a contour C that encloses the 
points uj = iujni n 7^ 0: 



duo 
c 2m 



log- 



-ioj, k) 



k'^ + 3m^ 



(44) 



where e^'^ (17 —> 0"*") is a convergence factor. The logarithm in the integrand can be expressed 
as the difference of two logarithms. The integral of the first logarithm reproduces the sum of 



the n 7^ terms in (^Bj). In the contour integral of the log(fc^ + Bm^) term, the contour can 



be deformed to wrap around the pole at = of the Bose-Einstein factor. By the residue 
theorem, this reproduces the n = term in (|43|) . 

The integrand in has a logarithmic branch cut in uj that runs from —ujL{k) to 
+ujL{k)^ where 00 ^{k) is the quasiparticle dispersion relation for longitudinal gluons. This 
dispersion relation satisfies k"^ — YlL^—iuJi, k) = Q 01 







k"^ + 3m^ 



lol , UJL + k 

1 - TTT log T 

2k lol — k 



(45) 



The integral also has a branch cut in u running from —k to +k due to the function 11^. 
This branch cut represents the effects of Landau damping. We choose both branch cuts to 
run along the real axis so that they overlap in the region —k<uj<k. 

The contour C can be deformed into a contour that wraps around the branch cuts from 
as in Fig. ^, and the contour can then be collapsed onto the branch cut. 

We identify the contributions from k <\uj\ < uiiik) as the quasiparticle part of J-'l and 
denote it by J^L,qp'- 



L,qp 



Tlog 



-I3k 



+ ^i^L- k) 



(46) 



We identify the contribution from < as the Landau-damping part of Tl and denote it 
by J^LM- 



^LM = - f f duj (pL 
H Jk Jo 



1 



1 



30^-1 2 



(47) 



The angle 0l is 



0i(ct;, k) = arctan 



2 '"-g k 



k"^ + 3m^ 



(48) 
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FIG. 4. The longitudinal free energy can be expressed as an integral around a contour C that 
wraps around the branch cuts of log(A:-^ — IIj^). 



where L = log [{k + uj)/{k — oj)]. The angle (pii^, k) vanishes as ^ and as — ^ oo and 
its value is vr at u; = /c. The denominator of the argument of the arctangent in (|4^) vanishes 
at some point uo = uj*{k). The arctangent in (^) is the principal branch for Q < uo < uj*{k), 
and it jumps to the next branch for u}*{k) < u < k, so that (pL remains continuous. 

The complete longitudinal free energy J-'l is the sum of the quasiparticle term ( ^Uf ) 
and the Landau-damping term (^71) . The quasiparticle integral is convergent, because ul 
approaches as a gaussian at large k. The Landau- damping term has ultraviolet divergences 
that are regularized by the (i-dimensional integral over k. In order to extract the divergences 
analytically, we make subtractions in the integrand that make the integral convergent in 
d = 3, and we then extract the poles in d — 3 from the subtracted integral. Our choice of 
the subtracted integral is 

2A;3 4A;(P + 3m2) ^8A;2(P + 3m2)2 J ■ ^ ^ 

After subtracting this from (^), we can take the limit d ^ 3: 

^ poo poo 

+ — / duj k'^. 

If we impose ultraviolet cutoffs < A and < A on the energy and momentum, the 
subtraction integral ( ^9] ) has a power divergence proportional to m^A^ and logarithmic 
divergences proportional to m^T^ log A and m^logA. The divergence proportional to m?gh? 
cancels against the corresponding divergence in the transverse Landau- damping subtraction 



n JkJo 



3Timgijj 



2k^ 



37!'mgLJ 



4, ,2 



2k{k^ + 3ml) 4A;2(A;2 + 3m2)2- 



(50) 
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integral pO|). Dimensional regularization throws away the quadratic divergence and replaces 
the logarithmic divergences by poles in e. In the limit d ^ 3, the individual integrals in (^91) 
are given by ( [A.7| )-( [AI9D in the appendix. The result is 



,7-(sub) 
L,Ld 



TT 



3 ,(3ml' 



1 C(-l) 

- + 2^^-21og(2vr) 



log2-l , ^ 1' 



(51) 



Our final result for the longitudinal free energy is the sum of 



m and (51). 



D. Renormalized Free Energy 



To complete our calculation of the HTL free energy, we must determine the counterterm 
AjF in (|l^) that cancels the divergences in {d — 1)J-'t + J^l- The ultraviolet divergences have 
been isolated in the subtraction terms (|39| ) and ( |4^ ) for J^t and (|5lD for J^l- The poles in 
e proportional to m^T^ cancel between J^t and J-'l. The cancellation can be traced to the 
fact that the logarithms cancels in the following combination of self-energy functions: 



7 O ~l~ 



-iuj^ k) 



3mt 



fc2 



k 



2 ■ 



(52) 



The remaining divergence is proportional to m^. In the minimal subtraction renormalization 
prescription, AjF is chosen to cancel only the pole in e and no additional finite terms: 

9 



1287r2 



-m„ 



(53) 



If we use a momentum cutoff A, there is also a quadratic divergence proportional to m^A^ 
coming from the transverse quasiparticle subtraction term (pTj). This would be canceled by 
an additional counterterm in ( ^3]) proportional to m^A^. 

To obtain the final result for the HTL free energy, we insert into (|18[) the sum 
of (H), (H), (0) and (H) for J^t, the sum of (HD, and for Tl, and (H) for 
AjF. We can simplify the result by evaluating the temperature independent integrals nu- 
merically. Since the only scale in these integrals is the mass rrig, they reduce to multiplied 
by a numerical coefficient. The contributions to from the second integral in ( |5^ ) and 
the second integral in (|4T|) are —8.656 x 10~^m^ and —4.554 x 10~^m^, respectively. The 
contribution to J-'l from the integral of the second term in (|46| ) and from the second inte- 



gral in (0) are 1.225 x 10' 
renormalized free energy is 



-2m^ and 4.290 x 10" 



m 



5' 



respectively. Our final result for the 



HTL 




21og(l - e" 



+ log 



1 - e"^^^ 
1 - e-/^^ 



„ , du —5 

27r3 7o e^'^ - 1 

H m„i H -m„ 

16 » 647r2 5 



A; (i/c \(hr, — 



log ^ -0.332837 



(54) 
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where /is = y' Aire'"' fi is the renormahzation scale associated with the modified minimal 
subtraction (MS) renormahzation prescription and 7 is Euler's constant. The dispersion 
relations uoT{k) and uJiik) are the solutions to ( |32D and (^S]). The angles (pTi'^-.k) and 
(j)L{uj,k) satisfy ( p5|) and 



IV. HIGH-TEMPERATURE EXPANSION 



If the temperature is much greater than the gluon mass parameter nig, the HTL free 
energy can be expanded in powers of rrig/T. Such an expansion is contrary to the spirit of 
HTL perturbation theory, which is an expansion in powers of g and 6 with fixed. It is 
nevertheless interesting, because it can be compared directly with the perturbative expansion 
for the QCD free energy. In this section, we compute the high-temperature expansion for 
the HTL free energy through order (m^/T)^. 



A. Separation of Scales 

The high-temperature expansion for the HTL free energy could be deduced directly 
from the final expression (^l)). However, this is difficult because of the delicate cancellations 
between quasiparticle and Landau-damping terms and between transverse and longitudinal 
contributions. It is easier to carry out the high-temperature expansion starting from dimen- 
sionally regularized contour integral expressions for J^t and !Fl- If we separate out the free 
energy of a massless boson given in (p6D , the transverse free energy (^) becomes 

Since the logarithm goes to zero as 00, there is no need for a convergence factor e'''^. 
We can choose the logarithmic branch cut in ( p5D to run from —uxik) to —k and from +k 
to +ujT{k). The contour C can then be deformed so that it wraps around the branch cuts 



that run along the real axis between —oj^ik) and +{jj'j'{k). In the expression (|43|) for the 
longitudinal free energy, the first term can be evaluated analytically using ( |A.10| ). In the 
second term, it is convenient to subtract log fc^ from the integrand before writing it as a 
contour integral as in (H5). The resulting expression is 



3^ ,!/■/■ '^^ i„„ A Hi(-iu;,A;)~ 



Tl = -—mlT + - / log 1 ^\ ' ' \ , (56) 

where the contour C wraps around the branch cuts that run along the real axis between 
—ujiik) and +ujL{k). 

The integrals in ( |55| ) and (^Bj) involve two energy and momentum scales: the "hard" 
scale T and the "soft" scale irig. The terms in the high-temperature expansion can receive 
contributions from both the hard scale and the soft scale. Dimensional regularization makes 
it easy to separate these contributions. We can obtain the soft contribution by expanding 
the Bose- Einstein factor in powers oi uj/T: 
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where i?„ are the Bernoulh numbers. The only scale in the resulting dimensionally regular- 
ized integral is rrig. Thus, upon inserting ( pT]) into (^) and (^) , we obtain expansions in 
powers of rrig/T. These are the soft contributions to the high-temperature expansion. We 



can obtain the hard contributions by expanding the logarithms in ( p^D and ( |5BD in powers 
of or, equivalently, in powers of 11^ and 11^: 




n 



-iuj, k) ' 



Uii-iuj, ky 



A;2 



OO 1 

E- 

OO 1 



n=l 



n 



n 



-iuj, k) 



ni(-ia;, A;) 



fc2 



(58) 
(59) 



The resulting dimensionally regularized integrals involve only the scale T. Thus, upon 



inserting (|58|) and (|59|) into (|55D and (|56D , we obtain expansions in powers of m^/T^. These 
are the hard contributions to the high-temperature expansion. In the next subsections, we 
calculate the expansions for the soft and the hard contributions through order (m^/T)^. 



B. Soft Contributions to Tt and 

The soft contributions to the high-temperature expansion of the free energy consist of 
the first term in (|56|) and the soft contributions from the second terms in (|55|) and (p6|). 
We first consider the soft contribution to the transverse free energy. After inserting the 
expansion into (|55D , there are no longer any poles on the imaginary uj axis. Having 
chosen the branch cuts to lie on the real axis in the interval —ujt < uj < ujt, the contour 
C can be taken to infinity in all directions. Making the change of variables z = l/u, the 
contour becomes a counterclockwise circle around z = 0. The soft contribution to J^t can 
then be written as 

(60) 

The logarithmic factor is an even function of z that is analytic at z = 0. It can therefore be 
expanded as a power series in z'^: 

iosfi-?FfeM) = EJ-.K*)^*. (61) 




{1/zy - p 



p=i 



where Tp{m, k) is a polynomial of degree p in and k'^. Inserting this expansion into (^ 
the contour integral vanishes unless n = 2p. The soft contribution (^) then reduces to 



-I OO 73 „ 

_ A V p^T'-^P / T (m, k) . (62) 
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However, since the integrand Tp{m, k) is a polynomial in k, the dimensionally regularized 
integral over k vanishes. The soft contribution to J^t is therefore zero. 

We next consider the soft contributions to the longitudinal free energy coming from the 
second term in (|56D . Performing the change of variables z = 1/lo yields 

where the contour is a counterclockwise circle around z = 0. The logarithmic factor is an 
even function of z that it is analytic at z = 0. It therefore can be expanded as a power series 

2 

m z . 

log (l - MZ^^ = £ L,{m, k)z'^ , (64) 

where Lp{m, k) is a polynomial of degree p in and k"^. The contour integral in ( |63D 
vanishes unless n = 2p. But then the integral over k vanishes in dimensional regularization 
because Lp is a polynomial in k. Thus the soft contribution to the second term of (^) is 
zero. The only contribution to the free energy from the soft momentum scale rrig is therefore 
the first term in which comes from the n = Matsubara mode. 



C. Hard Contribution to !Ft 



The hard contributions to the high-temperature expansion of J-'t consist of the term 
in (^) and the power series in obtained by inserting (^) into The first term in 

that expansion is 



^(1) 

•J J-" 



1 



duj IlT{—iuj, k) 1 



2 A Jc 2m uj^ - fc2 ef^^ -1 
The function multiplying the Bose-Einstein factor can be written as 



(65) 



IlT{—iuj,k) 3m 
a;^ — /c^ 



UJ 



2F 2k 



UJ UJ + k 



UJ — k 



(66) 



Thus the integrand in ( |65|) is the sum of two terms, one with poles aX uj = ±k and the 
other with a branch cut that runs from uj = —k to uj = k. The integral of the pole terms 
can be evaluated using the residue theorem. In the branch cut term, the contour can be 
collapsed onto the branch cut, and then expressed as an integral over positive values of uu. 
The resulting expression is 



1 



k \eP^ - 1 



+ 



I fk /I 1 

— / duj UJ \ — \ — 

A;3 70 \eP^ - 1 2 



(67) 



The integrals can be evaluated analytically. In the last term, the integrand is evaluated by 
integrating first over k and then over uj. The final result is 
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= I - 



TT 



, 2r7-i2 



2 



1 , C'(-i^ 



The second term in the expansion of Tt in powers of is 

2 



2 log(27r) - 2 



(68) 



^(2) 



1 r r duj IIlT{-iuj,k) 
4 A ic 2^ I _ p 



1 



(69) 



Using (|66D , we can separate the integrand from (^) into a double pole term, a double 
logarithm term, and terms that are products of a single pole and a logarithm. After using 
the residue theorem and collapsing the contour onto the branch cut, the expression (|69D 
reduces to 



^(2) 



9 4 /• r3 + 2e 
32 ^ 



1 

1 + 2 



Ve^'^ - 1 



2 ^ 3 ^ 

Jo duo Ve^^ 



1 



log 



k + u 



(70) 



The double integrals can be evaluated by integrating first over k and then over u. Expanding 
around e = and keeping terms only through order e°, the final result is 



(2) 



3 JT'' 



(2vr) 



log2-5 /I ^, ^ ^ 1\ 7 TT^ 



(71) 



D. Hard contribution to J^l 



The hard contributions to the high-temperature expansion of Tl can be obtained by 
inserting the expansion (BO) into (p§). The first term in that expansion is 



L 



1 



doj Yiii—i^, k) 1 



2 A Jc 2-111 



fc2 



e/3<^ _ 1 ' 



(72) 



Inside the contour C, the integrand has a pole at = from the Bose-Einstein factor and 
a branch cut in uj running from —k to k from the Hl factor. The residue of the pole is 
Since there is no scale in the subsequent integral over momentum, the integral 
vanishes using dimensional regularization. The pole can therefore be ignored. Collapsing 
the contour onto the branch cut, the expression reduces to 



(1) 



3 2 ^fi 

-m- — — 

2 ^{2tiY 



dk k 



f 

Jo 



duj UJ 



1 1 

T+2 



(73) 



Integrating first over k and then over u, our final result is 



TT 



2rri2 



fin 



(27r) 



1 , C'(-i' 



21og(27r) 



(74) 
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Multiplying the corresponding transverse term (|68D by 2 — 2e and adding to (^), we see 
that the 1 /e poles cancel: 

The second term in the expansion of JF^ in powers of is 



(75) 



^(2) 
■J L 



I r r duo f IlL{-iuj , k)' 
4 A ic 2^ V 



1 



(76) 



Inside the contour C, the integrand has a pole at = and a branch cut in u running 
from —k to k. Again the contribution from the pole vanishes after integrating over k using 
dimensional regularization. The contour can therefore be collapsed onto the branch cut 
and (|76|) reduces to 



-m„ / TT / duj uj 
2 ^kk^ Jo 



UJ ^ k + UJ 
2k ^ k — UJ 



1 



1 



(77) 



This can be evaluated by first integrating over k and then over uj. Expanding around e = 
and including terms through order e°, this reduces to 



-77(2) 
L 




ft. 



(27r)' 



;i-log2)('i-21og(27r) + 27-^) 



Multiplying the corresponding transverse term ([7T|) by 2 — 2e and adding to 



-(2) 



9 ,/r2- 

32 



(2vr)' 



3 + 27 - 2 log(27r) 



(78) 
we obtain 
(79) 



E. High-temperature expansion for ^htl 



The high-temperature expansion for the HTL free energy is obtained by inserting the 
high-temperature expansions for J-'t and J-'l into ([I8D along with the counterterm (0). The 
expansion includes the term from (|55D, the m^T term from (|56D, the m?T'^ term from ([75D, 



and the term from (|79D . Combining all the terms, the renormalized free energy is 



HTL 



ideal 



1 



3m2 ^ 
47r2T2 



30 



45 / 



' 3ml ' 
47r2T2 



1.232 



3/2 



/ 3m2 \ . fi , fl. 



(80) 



In the high-temperature limit, the thermal gluon mass parameter approaches = ^a^T^. 
If we make this identification, the HTL free energy is a selective resummation of terms in the 
QCD free energy that are higher orders in as- The expansion parameter Sml/An'^T'^ in (80) 



coincides with ag/Tr. Thus the high-temperature expansion ( pOD can be compared directly 
with the expansion of the QCD free energy in powers of a]^'^ which is given in (|l|). Comparing 
the coefficients of the expansions (|80|) and (|l]), we make the following observations: 
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The HTL free energy overincludes the correction by a factor of three. In conven- 
tional perturbation theory, the correction — ^a^/vr arises from interactions between 
massless transverse gluons. In HTL perturbation theory, this correction is separated 
into two terms: — ^as/vr from the masses of the quasiparticles, which is included at 
leading order, and -|-^as/vr from interactions between quasiparticles, which is included 
at next-to-leading order. An explicit calculation of the next-to-leading order correction 
is in progress [ffH . 



The a^/^ term, which arises from the effects of electric screening, is included exactly at 
leading order in HTL perturbation theory, provided that the thermal mass parameter 
is chosen to be = AnasT"^ up to corrections of order a^. 



• If we choose the renormalization scale /X3 in ( PU] ) to be of order irig, the HTL free 
energy includes a fraction ^ of the log as correction. 

In the HTL free energy, the oversubtracted as term combines with the terms of order a^ and 
higher to give an overall correction that is negative in spite of the large positive correction 
from the a^J'^ term. 



V. LOW-TEMPERATURE LIMIT 



In this section, we derive the low-temperature limit of the HTL free energy for fixed rUg. 
QCD undergoes a phase transition to a confining phase as the temperature decreases below 
some critical temperature Tc. Because of this phase transition, we do not expect J-htl in 
the limit T <^ nig to bear any resemblance to the free energy of QCD at T < Tc. However, 
if J-'nTL is to be used as a phenomenological model for J^qcb for T > Tc, it is worthwhile to 
determine the qualitative behavior of JFhtl for extreme values of its parameters. 

In the low-temperature limit, J-'utl is proportional to m^. The coefficient could be 
extracted directly from our final expression (|5^, but it is simpler to compute it directly 
from our original expression for the transverse and longitudinal free energies. 



A. Transverse Free Energy 



Our original expression (^) for the transverse free energy involves a sum over the discrete 
Matsubara frequencies Un = 27mT. As T 0, the sum approaches an integral over the 
continuous euclidean energy u: 



1 

47r 



duj / log F-Fa;^-FHT(cu,A;) 

00 '- 



(81) 



Since Hy is a function of the combination ujjk only, it is convenient to rescale the energy 
variable \>y uj kuj. Integrating over the angles of k and using the fact that the integrand 
is an even function of cj, the integral reduces to 



— ^^Ai / duj 



dkkHog (1 + cu^)^:^ + Ht(w, 1) 



(82) 
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The dimensionally regularized integral over k can be evaluated analytically usin^ 
giving 



1 n^n'-Y^) „s-. 



27r(27r)'^ d + 1 
Expanding around d = 3, we get 



-/i 



d+l 



(83) 



327r ^ V V / (2vr)'=' I Ve ' 2 



1 1 \ /""^ 

du - / rft^ f^iuj) log /T(t^) 

JO 



where the function in the integrand is 



arctan u 



1+UJ 



2 ■ 



The first integral in (p^ ) can be evaluated analytically, but the second integral 
evaluated numerically. The final result is 



_9_ ' 



(in. 1,^0.200871 



(85) 

must be 

(86) 



B. Longitudinal Free Energy 



In our original expression (^) for the longitudinal free energy, the sum over the Mat- 
subara frequencies Un approaches an integral over u in the limit T ^ 0. Integrating over 
the angles of k and rescaling the energy variable by — kuj, the expression reduces to 



3-d r 

2-n {27iY^ Jo 



duj / dk riog r -nL(tu,l) 

Jo L 



(87) 



The integral over k can be evaluated analytically using ( |A.10| ). Expanding around d = 3, 
our expression reduces to 



9 



'3ml 



1 1 
7+2 



Jo 



{2nY 

where the function in the integrand is 

= 1 + u! ^arctan — — ^ . (89) 

Evaluating the first integral in (|88|) analytically and the second numerically, we obtain: 

^d 



9 J3m'; 
m„ 

StT M 



(2vr)^ 



7r(l-log2) ^11 
2 



0.176945 



(90) 
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C. Renormalized Free Energy 



We can obtain the T — limit of the HTL free energy by inserting (|8^) , (0), and the 
counterterm into the expression (IIS) for J-'utl- The result is 



9 / m„ \ 
_ log _1_ 0.332837 ™ 



(91) 



which is identical to the term in our final expression (p4|) for the renormalized free energy. 
This is an efficient way of computing the constant under the logarithm in that term. 

Among the next-to-leading order terms in the low-temperature expansion for J-rtl is the 
m^T^ term in (0). To identify the remaining next-to-leading order terms, we consider the 
T — >■ limit of each of the integrals in (0). The first integal in (0) is the sum of a transverse 
term involving ujt, a longitudinal term involving ul, and a subtraction term. The integral 
of the subtraction term gives a contribution proportional to T^, and thus contributes only at 
next-to- leading order in T^/m^. The transverse and longitudinal terms are dominated by the 
region of small momentum k. The dispersion relations in this region can be approximated 
by 

Mk) ^m, + (92) 

Mk) ^ ™, + 1^^. (93) 

Making the approximation log(l — e~^'^) ~ —e~^'^ in the integrand, we see that both the 
transverse and the longitudinal terms reduce to Gaussian integrals in k. Their contributions 
to the free energy scale as T {nigT)^^'^ ^-'^a/T ^ which falls faster than any power of T. 

The Landau-damping term in (|5^) involves the angles l and 0t defined in (^) and (^8|) . 
The Bose-Einstein factor in ( ^4] ) constrains uo to be of order T. In the integral over k, there is 
a logarithmic ultraviolet divergence proportional to that cancels between (p^ and 0t- The 
corresponding infrared cutoff is of order nig for the 01, integral and of order {nigUjY^'^ for the 
0T integral. Thus the leading contribution from the Landau-damping term is proportional 
to m^T^ log(mg/T). 

The low-temperature limit of J-rtl is sensitive to the value of /is. In particular, the 
coefficient of the term in ( |9TD changes sign at ~ 0.117 nig. For larger values of /is, 
the ratio J-'-rtl / ^ideai approaches +oo hke as T — > with nig fixed. For smaller 

values of /is, the ratio approaches — oo like For /is = /is, the divergence is less severe 

approaching — oo like 



VI. THERMODYNAMICS 

In this section, we examine the thermodynamic functions at leading order in HTL pertur- 
bation theory. We first derive a convenient expression for the trace anomaly density £ — 3V 
in terms of derivatives of the parameters nig and /is with respect to T. We give a prescrip- 
tion for the temperature dependence of nig which is motivated by the high-temperature limit 
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of QCD. We then compare our leading order HTL calculations for the pressure V and for 
£ — 3V with results from lattice simulations of pure-glue QCD and with gluon quasiparticle 
models. 



A. Energy Density 



Once the free energy JF(T) is given as a function of T, all other thermodynamic functions 
are determined. In particular, the pressure V and the energy density £ are 



dT 
IT' 



The combination £ — 3V can be written as 



£ -3V = -T 



d (T 



dT VT4 



(94) 
(95) 

(96) 



This combination is proportional to the trace of the energy-momentum tensor. In QCD with 
massless quarks, it is nonzero only because scale invariance is broken by renormalization 
effects. We will call it the trace anomaly density. It of course vanishes for an ideal gas of 
massless particles. However, it also vanishes for a gas of quasiparticles whose masses are 
linear in T and whose interaction are governed by a dimensionless coupling constant that is 
independent of T. 

Our expression (|5^) for the HTL free energy is a function of three variables: T, m^, and 
/is. Thus the temperature dependence of Twvl is determined only after the T-dependence 



of mg{T) and iJi^{T) is specified. Using the chain rule, the expression (^) for £ — 3V can 
be written as 



£'KTL — 3Phtl 



d dirta d dii-i d 



dT dT dnig dT dfxs 



HTL 



(97) 



But J^htl/^^ is a homogeneous function of degree zero in the variables T, m^, and /X3, and 
it therefore satisfies 



d d 



^1 



HTL 



2^4 



0. 



This can be used to eliminate the partial derivative with respect to T from (| 

d /, m„\ d ^ „ d 



(98) 



£}1TL — 3V] 



HTL 



_T^(log!^ 



"^^d^g^^^^ - ^dT f ) ^W/^^"^- 



(99) 



This expression vanishes if rrig and /is are exactly linear in T. The partial derivative with 
respect to /is in ( pUj ) is simple. The partial derivative with respect to rrig is more complicated 
because in addition to the exphcit dependence on nig, there is the implicit dependence of 
then functions Ut, ujl, (pT, and (pL on m„. 



26 



B. T-dependence of nig and ^3 



Our leading order results for the thermodynamic functions depend on the thermal gluon 
mass parameter rUg and on a renormalization scale fi^. These parameters are completely 
arbitrary in the sense that the dependence on them will be systematically canceled by higher 
orders in HTL perturbation theory. Since they are being used as a device to reorganize the 
perturbative series for thermal quantities, they should depend on the temperature T. If 
higher order calculations in HTL perturbation theory were available, reasonable values of 
mg(T) and fJ^siT) could be determined by the condition that the higher order corrections be 
well behaved. In the absence of such calculations, the best we can do is to make physically 
motivated estimates of mg{T) and n^iT). 

A useful source of intuition on reasonable values of rUg is the high-temperature limit of 
QCD. In this limit, the expression for the thermal gluon mass is 



m%T) = ^aMT', (100) 



where ^4 is a renormalization scale proportional to T. If we are to use (|10CI|) , we must 
also specify the T-dependence of n^. The expression ( |100| ) can be derived in the form of a 
sum-integral over euclidean four- momentum 0: 

>nl{T) = '-^aAp,)f^'^, (101) 

where = /cq + k^. The only momentum scales in the sum-integral are integer multiples of 
the lowest Matsubara frequency 2iiT. This suggests that a reasonable value is n^iT) 2'kT. 

If we use a parametrization of the running coupling constant as{ii4) that diverges at 
some small momentum scale, our prescription (|100|) for the thermal gluon mass will diverge 



at a sufficiently small temperature. A parameterization of the running coupling constant 
that includes the effects of two-loop running is 

^^U- , ,102) 



where L = log (^jjl/Ay^j. This expression ( |102| ) diverges at ^4 = Aj^, so the thermal gluon 



mass rrigiT) in ( p.OCI|) diverges at the temperature T given by n^iT) = A^. 

In their original paper on screened perturbation theory [jl2|, Karsch, Patkos, and Pe- 
treczky chose the mass of the scalar quasiparticle to be the solution of a one-loop gap 
equation. An analogous choice for the gluon mass parameter would be to take mg{T) to be 
the solution of the integral equation 

" dk k ^^^^^_^ . (103) 



In the limit Q;s(/i4) -C 1, the solution reduces to ( |100|) , with the leading correction propor- 
tional to a^/^. This correction would contribute an additional a^/^ term to J-'utl through 
the term in (pOt). It would be canceled at next-to- leading order in HTL perturbation 
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theory by a similar contribution from the term from the one-loop diagram with an HTL 
counterterm. In order to get the correct a^^^ term at leading order in HTL perturbation 
theory, must be given by (|100|) up to corrections of order af. We choose here to use the 
simplest possibility ( |100D . 

We next consider the renormalization scale fi^. The dependence on fi^ arises because the 
one- loop HTL free energy has a logarithmic ultraviolet divergence proportional to rrig. The 
divergence has a three-dimensional origin: it arises as the momentum k ^ oo with fixed 
energy u. Similar divergences proportional to g'^T'^rrig and g'^T'^ will appear at next-to- 
leading order and next-to-next-to-leading order in HTL perturbation theory. The infrared 
cutoffs on the logarithmically divergent integrals are provided either by the energy u, which is 
of order T, or by the thermal gluon mass nig. If we use the weak-coupling limit ( |10CI| ) for m^, 
the divergences proportional to m^, g'^T'^m'^, and g'^T^ will cancel exactly, leaving logarithms 
of T/rUg. These logarithms give the a^loga^ terms in the weak-coupling expansion (|I]) for 
the QCD free energy. 

The high-temperature limit of the HTL free energy is rather insensitive to the value of 
yU3, but the low-temperature behavior is very sensitive. With the parameterization (|102|) 
for the running coupling constant, the thermal gluon mass diverges as T approaches the 
value where fXiiT) = A-^. The HTL free energy can therefore be approximated by the low- 
temperature expansion given in (pi]). If we choose /is to be proportional to irig, J-'nTh/^ideai 
diverges like m^(T), approaching +00 for /X3 > fi^ and —00 for fj,^ < /Xg, where /ig = 
0.717mg. If we choose /is = /tg, the HTL free energy diverges more slowly as rn?g{T). We 
can minimize the pathological behavior of J^htl at low temperatures by choosing the value 
/is(T) ^ Q.inmgiT). 

C. Comparison with Lattice Gauge Theory 

Lattice gauge theory can be used to calculate the thermodynamic functions of QCD 
from first principles, with all errors under control. We can therefore compare our leading 
order HTL results with the correct answer provided by lattice calculations. Boyd et al. |^ 
have calculated the thermodynamic functions of a pure SU (3) gauge theory for temperatures 
up to about 5Tc. They used Monte Carlo simulations with high statistics on a lattice as 
large as 32^ x 16 to calculate plaquette expectation values for about 20 points in T ranging 
from to about 5Tc. They extracted the trace anomaly density £{T) — 3V{T) directly from 
the plaquette expectation values at and T. The pressure V(T) was extracted from the 
plaquette expectation values at temperatures between and T by computing an integral. 
Their final results were obtained by extrapolating to the continuum limit and were presented 
in the form of continuous interpolation curves. The results of Boyd et al. for V{T) and 
S{T) — 3V(T) are shown in Figs. || and ^ The pressure is normalized to the pressure "Pideai 
of an ideal gas of massless gluons. The trace anomaly density is normalized to the energy 
density £^ideai = SPideai of the ideal gas. The normalized energy density S/Sideai is given by 
the sum of the values in Figs. ^ and ^ The interpolation curves are represented by diamonds 
at a set of discrete points. The size of the diamonds indicates the typical error of about 3% 
at any single value of T. 

The leading order HTL results for the pressure and the trace anomaly density are shown 
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FIG. 5. The pressure for pure-glue QCD as a function of T /T^. The HTL free energy is shown 
as a shaded band that corresponds to varying /i3 and /X4 by a factor of 2 around their central 



values. The weak-coupling expansions through order as 



a 



3/2 



a^, and a 



5/2 



are shown as dashed 



lines labeled by 2, 3, 4 and 5. The lattice results of Boyd et al. H are shown as diamonds. 



in Figs. ^ and |^. We use the expression ( |100|) for the thermal gluon mass, with the running 
coupling constant given by ( |1U2[ ) and with Tc = 1.03Aj;jg. To illustrate the sensitivity to the 
choices of renormalization scales /^a and yU4, we take their central values to be 

/i3(T) = 0.717mg(T), (104) 
/i4 = 27rT, (105) 

and we allow variations by a factor of 2 in the coefficients on the right side of (|104|) and ([1051) . 
The shaded bands in Figs. ^ and ^ indicate the resulting range in the predictions. The 
ranges are dominated by the variation in /i4 at the largest values of T shown and by the 
variation in /is for T near T^.. For comparison, the predictions from the QCD weak-coupling 
expansion (|Tp with /i4 = 27rT are also shown in Fig. H and Fig. ^. The expansion (Q) of 
the pressure truncated after orders a^, a^''^; ^^"^ shown in Fig. || as the dashed 

lines labeled 2, 3, 4, and 5. The corresponding predictions for the trace anomaly density are 
shown as dashed lines in Fig. ^ They were obtained by evaluating the derivative in (|96D 
numerically. 

The HTL pressure shown in Fig. |] has the correct shape at the largest values of T, but 
the slope also remains small for T near Tc. The small slope can be understood from the 
fact that our expression for mg{T) in ( |100| ) is almost linear in T. If it were exactly linear in 
T, then Phtl/^^ would be constant. The factor of the running coupling constant in (|100|) 
causes mg(T) to grow a little slower than Unear in T, which causes Vutl/T^ to increase 
slowly with T. 

The leading order HTL free energy in Fig. ^ lies below the lattice results at the highest 
temperature available. This deviation has the correct sign and roughly the correct magnitude 
for the inclusion of the next-to-leading order correction in HTL perturbation theory to give 
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1 2 3 4 5 

T/To 

FIG. 6. The trace anomaly density for pure-glue QCD as a function of T/Tc- The HTL trace 
anomaly density is shown as a shaded band that correspond to varying /is and fi^ by a factor of two 
around their central values. The result from differentiating the weak-coupling expansions through 
order Ug, o^J"^ , a^, and aV^ are shown as dashed lines labeled 2, 3, 4, and 5. The lattice results 
ofBoyd et al. Q are shown as diamonds. 

a better approximation to the pressure. The next-to-leading order correction to VHTifPideai 
should be positive at large T, since it must approach +Yas/7r at asymptotically large 
temperatures. 

The HTL prediction for the trace anomaly in Fig. ^ is very small. The reason for this is 
that the differentiation in (^) increases the size of the a^^^ term relative to the term by 
a factor of 3/2. As a consequence there is a near cancellation between the over subtracted 
term and the a^/^ term. If the next-to-leading order correction in HTL perturbation theory 
is dominated by the remaining as term, it would give a negative contribution, increasing the 
discrepancy with the lattice results. However, the fact that there is a near cancellation at 
leading order makes it dangerous to predict the effect of the next-to-leading order correction. 



D. Comparison with Quasiparticle Models 



The HTL free energy ( |54D can be interpreted as essentially that of an ideal gas of gluon 
quasiparticles, except that it is modified in such a way as to consistently take into account 
the screening effects of the plasma. It is worthwhile to compare it with the phenomenologi- 
cal quasiparticle models that have been used to describe QCD thermodynamics [|T^,|r^. In 
the most recent analyses by Peshier et al. [T^ and by Levai and Heinz |T^, the quasipar- 



ticle model is an ideal gas of transverse gluons with a temperature-dependent mass. Their 
expressions for the pressure and energy density are 



^ = i^^/„ * * 7Fhi:^-^<^'- (106) 
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S = ^l dkeur^^^^+B{T). (107) 

The "bag function" B{T), which cancels in the entropy density S = {S+V) /T, is determined 
up to an integration constant Bq by the thermodynamic consistency condition 

dV 

S = T—-V. (108) 
The transverse dispersion relation ujt was taken to be 

4 = + Ik (109) 
which is the solution to (B^) in the limit k ^ 



"5- 

We now compare the expression ( |106| ) for the pressure in the quasiparticle model with 
the HTL pressure which is the negative of (p^ . The transverse quasiparticle contributions 
agree after an integration by parts. The only difference is the use of the simplified dispersion 
relation (|109|) instead of the solution to (|32D , but that is a good approximation in the 
most important momentum regions. In the quasiparticle model, the contributions from 
the plasmon, from Landau damping, and from zero-point energies are replaced by the bag 
function B{T), which is assumed to cancel in the combination S{T)+V{T). This assumption 
is not justified within the HTL framework. 

In their analysis, Levai and Heinz determined the mass mg{T) by fitting the lattice 
results for the pressure V{T) [|19[. Their result for mg{T) is approximately linear in T 
between 2Tc and 4.5Tc, with a coefficient of proportionality of about 0.75, and it rises to 
about 3.5T as T approaches T^. This is qualitatively similar to the behavior predicted by 
the expression (|100|) for the thermal gluon mass in the weak-coupling limit. If the scale 
/i4 in ( |100| ) is determined by fitting the lattice results at the highest available temperature 
near 4.5Tc, the result is fi^ ~ 18T. This is uncomfortably large for the scale of the running 
coupling constant. This large value of /i4 has a simple explanation from the point of view 
of HTL perurbation theory. The next-to-leading order contribution to V/Videai in HTL 
perturbation theory will be positive at large T since it must approach -|-Y«s/7r in the high- 
temperature limit. Fitting the lattice results without allowing for the next-to- leading order 
correction would require a smaller value of rrig and hence a larger value of n^. Thus the large 
value of /i4 obtained in quasiparticle models is an indication that the interaction between 
the quasiparticles cannot be neglected. 



VII. OUTLOOK 

We have calculated the free energy of a gluon plasma to leading order in HTL per- 
turbation theory. Extending the calculation to include quarks is straightforward. A more 
challenging problem will be to extend the calculations to next-to-leading order. This requires 
calculating two-loop diagrams in HTL perturbation theory. Such a calculation is essential 
in order to demonstrate that HTL perturbation theory avoids the convergence problem that 
plagues the conventional perturbative expansion. 
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In the high-temperature hmit, the two- loop HTL free energy will agree with the per- 
turbative expansion (|l]) up to errors of order of a^. The three-loop HTL free energy would 
agree up to order log as- In general, the n-loop calculation in HTL perturbation theory 
will include all the n-loop contributions from the scale T, which scale like {g'^)"'~^, and all 
the n-loop contributions from the scale gT, which scale like (7^^". One limitation of HTL 
perturbation theory is that it can not be used to calculate the contribution from magneto- 
static gluons with momenta of order g'^T, which first contribute to the free energy at order 
g^. These effects are inherently nonperturbative. It is possible that they could be calculated 
by a strategy analogous to that used by Kajantie et al. to compute the Debye mass for 
QCD. The magnetostatic gluons can be described by an effective theory, called magneto- 
static QCD (MQCD), which is a pure SU{Nc) gauge theory in three euclidean dimensions. 
The parameters of MQCD could perhaps be calculated by HTL perturbation theory, and 
then the effects of magnetostatic gluons could be calculated nonperturbatively by applying 
lattice gauge theory methods to MQCD. 

In HTL perturbation theory, the leading order quasiparticle dispersion relations are built 
into the propagator. At next-to-leading order in g, the quasiparticle dispersion relations 
and ujI have logarithmic infrared divergences proportional to g'^Tm that arise from magne- 
tostatic gluons. It would be unwise to build the next-to-leading order dispersion relations 
into the propagator, because this would introduce a sensitivity to the magnetostatic glu- 
ons into the one-loop free energy that would have to be canceled by higher-loop diagrams. 
The infrared divergences in the next-to-leading order quasiparticle dispersion relations will 
of course give a divergent contribution to the next-to-leading order free energy through 
two-loop diagrams that include a one-loop gluon self-energy correction. However, we ex- 
pect this divergence to be canceled by a divergence in the corresponding Landau-damping 
contribution. 

One of the advantages of HTL perturbation theory is that it can be applied to the real- 
time processes that are the most promising signatures for the quark-gluon plasma. With 
the exception of the production of hard dileptons , these signatures have been calculated 
only at leading order in ordinary QCD perturbation theory. There are two reasons to be 
concerned about the reliability of these calculations. One is that the higher order corrections 
for the signatures are probably at least as large and unstable as the higher order corrections 
for the free energy. There is therefore no way to determine the accuracy of the leading order 
calculation. The other reason for concern is that the conventional weak-coupling expansion 
does not give a good approximation to the equation of state. The equation of state is 
needed to infer a temperature T from the energy density of hot hadronic matter, and T 
is then used as a parameter in the calculation of the signatures. If the method used to 
calculate the signatures does not reproduce the equation of state, the whole framework is 
inconsistent. If the next-to- leading order calculation in HTL perturbation theory gives a 
good approximation to the equation of state, we can have confidence in the predictions for 
the signatures that are calculated within the same framework. 

HTL perturbation theory provides some justification for quasiparticle models of the 
quark-gluon plasma. However, it goes far beyond those models, because it provides a frame- 
work for systematically calculating the effects of interactions between the quasiparticles. 
This could have an enormous impact on the phenomenology of the quark-gluon plasma, 
because the physical picture that is suggested by HTL perturbation theory is dramatically 



32 



different from the conventional picture of the quark-gluon plasma as an almost ideal gas 
of ultrarelativistic quarks and gluons. In HTL perturbation theory, the quarks and gluons 
have thermal masses that are comparable to the temperature T, so they are only mildly 
relativistic. This dramatic difference cannot help but have significant phenomenological im- 
plications. There have been two previous studies of the effects of quasiparticle masses on 
the signatures of a quark-gluon plasma. Biro, Levai and Miiller have studied the effect 
of a massive gluon on the strangeness production in a quark-gluon plasma. They took the 
u and d quarks to be massless, they neglected the longitudinal mode of the gluon, and they 
took the transverse mode to have a temperature-independent mass of 500 MeV. They found 
that the gluon mass enhances the production of ss pairs at temperatures below 300 MeV. 
Additionally, the effects of thermal masses for quarks and gluons on charm production from 
a quark-gluon plasma has been studied by Levai and Vogt [Q. They neglected the longi- 
tudinal mode of the gluon and took the transverse gluons to have the dispersion relation 
a;|n = k"^ + m?g. They found that the thermal masses significantly enhanced the thermal 
production of charm at RHIC and at LHC 

The previous calculations of the effects of quasiparticle masses have serious theoretical 
inconsistencies, because introducing a transverse gluon mass by hand destroys gauge invari- 
ance. HTL perturbation theory solves this problem by introducing the transverse gluon mass 
in a gauge-invariant way. But the transverse gluon mass is intricately linked to other essen- 
tial aspects of relativistic plasma physics, including longitudinal gluons. Landau damping, 
and the screening of interactions. All of these effects are incorporated consistently within 
HTL perturbation theory. Thus it provides a foundation for developing a new phenomenol- 
ogy of the quark-gluon plasma in which the many-body aspects of the system play a central 
role. 
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APPENDIX A: INTEGRALS 



In this appendix, we collect the results for the integrals that are required to calculate the 
one-loop HTL free energy, the first few terms in the high-temperature expansion, and the 
low-temperature limit. We use dimensional regularization, so that ultraviolet divergences 
appear as poles in e. 

In the HTL free energy, the ultraviolet divergences are isolated in subtraction terms that 



must be expanded around e = 
quasiparticle subtractions are 



through order e . The integrals required to evaluate the 
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The corresponding integrals required to evaluate the Landau-damping subtractions are 
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(A.3) 
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There is also a T-dependent integral that multiplies a pole in e and must therefore be 
expanded to order e: 



duj 



UJ 
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(A.9) 



There are several nontrivial integrals required to compute the high-temperature expan- 
sion and the low-temperature limit for the HTL free energy. One of these integrals is 



dk log(fc2 ^ ^2N ^ M 2 2 

a + 1 



(A.IO) 



This integral with a = 2 is needed to compute the soft contribution to the free energy from 
the n = Matsubara mode. The same integral with a = 3 — 2e is needed to compute 
the low-temperature limit of the HTL free energy. In computing the (m^/T)^ term in the 
high-temperature expansion, we need the following infrared-divergent integral: 



°° , UJ 
duj 



— ^ = -(^- + 27-21og(27r; 



(A.ll) 



We also need the following integral which comes from changing variables in a momentum 
integral x = k/uj: 



dx x' 



-4-2e 



log 



a; + l_l + 21og2 7r2-81og2-10 
x-l ~ 3 ^ 18 ' 



(A.12) 



Since it multiplies an infrared pole in e from the integral ( [A.ll| ), it must therefore be ex- 
panded to order e. 
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